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The propagation of charm and bottom quarks through an ellipsoidal domain of quark gluon plasma 
has been studied within the ambit of non-equilibrium statistical mechanics. Energy dissipation of 
heavy quarks by both radiative and collisional processes are taken in to account. The experimental 
data on the elliptic flow of the non-photonic electrons resulting from the semi-leptonic decays of 
hadrons containing heavy flavours has been reproduced with the same formalism that has been used 
earlier to reproduce the nuclear suppression factors. The elliptic flow of the non-photonic electron 
from heavy meson decays produced in nuclear collisions at LHC and low energy RH1C run have also 
been predicted. 
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The nuclear collisions at Relativistic Heavy 
Ion Collider (RHIC) and the Large Hadron Col- 
lider(LHC) energies are aimed at creating a phase 
where the bulk properties of the matter are governed 
by (light) quarks and gluons. Such a phase of matter 
is called Quark Gluon Plasma (QGP). The study of 
the properties of QGP is a field of high current in- 
terest. Experimental results already available from 
RHIC li] have been successfully described by the ap- 
plication of relativistic hydrodynamics @, Q indicat- 
ing substantial collective dynamics in the system. In 
a hydrodynamic model one of the key assumption is 
that the system achieves thermalization. The ellip- 
tic flow parameter v^ F , is one such quantity which 
is sensitive to the early thermalization of the system. 

As the relaxation time for heavy quarks are larger 
than the corresponding quantities for light par- 
tons [4] the light quarks and the gluons get ther- 
malized faster than the heavy quarks. Therefore, 
the propagation of heavy quarks through a QGP 
(mainly containing light quarks and gluons) may 
be treated as the interactions between equilibrium 
and non-equilibrium degrees of freedom and the 
Fokker-Planck (FP) equation provides an appropri- 
ate framework |4Tjl3j for such studies. Since heavy 
quarks remain out of equilibrium i.e they are not a 
part of the equilibrated system and their production 
is restricted to the primordial stages of the collision, 
they can not decide the bulk properties of the sys- 
tem, rather act as an efficient probe to extract infor- 
mation about the system. Therefore, in the present 
work we will use the elliptic flow of heavy quarks [T^ — 
HU as a probe to extract the properties of QGP. 

In our earlier works (l7j . we have used the FP 
equation to study the experimental results from 
STAR and PHENIX collaborations on nuclear sup- 
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pression factor, Raa of the non-photonic electrons 
resulting from the semi-le pto nic decays of hadrons 
containing heavy flavours |l8l . Ii~9j . In the same ap- 
proach the Raa's have been predicted f20L [2l1 | for 
LHC and low energy RHIC runs [H [U In the 
present work we intend to study the elliptic flow 
measured at RHIC energy with the same formalism 
and input parameters that reproduce the Raa ■ The 
v^ F have also been predicted for low energy RHIC 
and LHC energy. 

The FP equation and the procedure adopted in 
evaluating the drag and diffusion co-efficients are 
discussed in detail earlier [T3, H3, Therefore, 
in the present work we outline the main results and 
refer to the earlier works for details. The evolution 
of heavy quarks momentum distribution function, 
while propagating through the QGP are assumed to 
be governed by the FP equation which reads, 
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where the kernels Ai and Bij are given by 



Ai — J d 3 kuj(p, k)k t 
Bij = J d 3 ku)(p 7 k)kikj. 
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for | p |— >■ 0, Ai — > jpi and By — > DSij where 7 and 
D stand for drag and diffusion co-efficients respec- 
tively. The temperature and the baryonic chemi- 
cal potential dependence of the transport coefficients 
enter through the thermal parton distributions (20j 
which appear in oj(p, k) as 



uj(p, k) 
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The basic inputs required for solving the FP equa- 
tion are the dissipation co-efficients and initial mo- 
mentum distributions of the heavy quarks. The drag 
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and diffusion coefficients have been evaluated by tak- 
ing in to account both the collisional and radiative 



spatial anisotropy as discussed above is given by: 



processes 
cone (H 



21]. In the radiative process the dead 



and Landau-Pomerenchuk-Migdal (LPM) 
effects are included, the details have been discussed 
in Ref. (2l|. In evaluating the drag and diffusion 



co-efficients we have used temperature dependent 
strong coupling, a s from [25]. The Debye mass, 
~ g(T)T is also a temperature dependent quantity 
which is used as a cut-off to shield the infrared di- 
vergences arising due to the exchange of massless 
gluons. 

The production of charm and bottom quarks in 
hadronic collisions is studied extensively 26]. The 
initial momentum distribution of heavy quarks in 
pp collisions have been taken from the NLO MNR 
code [27|. The results from the code may be tested 
by measuring the production cross sections of heavy 
mesons (containing c and b quarks) in pp collisions 
at RHIC and LHC energies. 

How can the heavy quarks probe the asymmetry 
of the thermalized system formed by light quarks 
and gluons? Consider a thermalized ellipsoidal spa- 
tial domain of QGP with major and minor axes of 
lengths l y and l x (determined by the collision ge- 
ometry) respectively. Now let us assume that a 
heavy quark propagates along the major axis then 
the number of interactions it encounters or in other 
words the amount of energy it dissipates or the 
amount of momentum degradation that takes place 
is different from when it propagates along the minor 
axis, because l y ^ l x . Therefore, the momentum 
distribution of electrons originating from the decays 
of charmed hadrons {D mesons) produced from the 
charm quark fragmentation will be anisotropic, thus 
reflecting the spatial anisotropy due to non-central 
collisions in the momentum distributions. The de- 
gree of momentum anisotropy will depend on both 
the spatial anisotropy and more importantly on the 
coupling strength of the interactions between the 
heavy quark and the QGP. In an extreme case of zero 
interaction strength the momentum anisotropy will 
be zero even if the spatial anisotropy is large. For a 
quantitative understanding we consider the follow- 
ing. The time evolution of the ith component of the 
average momentum for an ellipsoidal geometry can 
be obtained from the expressions: 



(Pi) ~ exp 



li/v 



7(i) dt 
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where v is the velocity of the heavy quark and 
li is the length along the ith direction. The ec- 
centricity of the momentum-space ellipse due to 
anisotropic momentum distribution resulting from 
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where e p is the eccentricity in the momentum space. 
The values of l y and l x are fixed by geometry of 
the collisions and can be estimated using Glauber 
model for different centralities. Both l x ^ l y and 
7^0 (non-zero interaction) are required for non- 
zero elliptic flow. The second Fourier coefficient of 
the momentum distributions is called elliptic flow, 
and we evaluate this quantity in the present 

work. 

The FP equation has been solved for the heavy 
quarks by using the Greens function technique 
(see [8| for details). The initial pr distribution is 
taken from NLO MNR formalism [27J ■ We convolute 
the solution (at the transition point i.e. when the 
QGP reverts back to hadronic phase) with the Peter- 
son fragmentation function [28[ of the heavy quarks 
to obtain the px distribution of the heavy {B and D) 
mesons. The momentum distribution, dN / dydpxd(j) 
of the non-photonic single electron spectra originat- 
ing from the decays of heavy flavour mesons - e.g. 
D —> Xev at mid-rapidity (y = 0) is obtained by fol- 
lowing the procedure described in [T3, [H, H3] ■ The 
coefficient of elliptic flow, v^ F then can be obtained 
as: 
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The system formed in nuclear collisions at rela- 
tivistic energies evolves dynamically from the initial 
to the final state. The time evolution of this systems 
has been studied by using the (2+1) dimensional hy- 
drodynamical model }3l| with boost invariance along 
the longitudinal direction [32j]. It is expected that 
the central rapidity region of the system formed af- 
ter nuclear collisions at LHC energy is almost net 
baryon free. Therefore, the equation governing the 
conservation of net baryon number need not be con- 
sidered for RHIC and LHC. 

The total amount of energy dissipated by a heavy 
quark in the QGP depends on the path length it 
traverses. Each parton traverse different path length 
which depends on the geometry of the system and on 
the point where its is created. The probability that 
a parton is produced at a point (r, cj>') in a QGP of 
ellipsoidal shape depends on the number of binary 
collisions at that point which can be taken as: 
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• PHENIX (0-20%) 



where i? is the nuclear radius and e is the eccentricity 
of the ellipse. Note that for central collisions (e = 0) 
Eq. [7] reduces to the expression for P(r, ft) used in 
Ref. Q for spherical geometry. A parton created at 
(r, ft) in the transverse plane propagate a distance 
L = \J R 2 — r 2 sin 2 (f)' — rcosft in the medium. We 
use the following equation for the geometric average 
of the integral which appear in the solution of the 
FP equation involving drag coefficient: 



/ rdrdftP(r, ft) J L/v cItj(t) 
/ rdrdftPir, ft) 



r = 



(9) 



where v is the velocity of the propagating partons. 
Similar averaging has been performed for the diffu- 
sion co-efficient. For a static system the temperature 
dependence of the drag and diffusion co-efficients 
of the heavy quarks enter via the thermal distri- 
butions of light quarks and gluons through which 
it is propagating. However, in the present scenario 
the variation of temperature with time is governed 
by the equation of state or the velocity of sound 
of the thermalized system undergoing hydrodynamic 
expansion. In such a scenario the quantities like T 
(Eq. [9]) and hence V2 becomes sensitive to velocity 
of sound (c s ) in the medium. 
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FIG. 1: Variation of v" F with p T at the highest RHIC 
energy. Experimental data is taken from [331 ]. 



The pt dependence of v^ F is sensitive to the na- 
ture of the initial distributions of heavy quarks, life 
time of the system, velocity of sound as well as the 
drag and diffusion coefficients of heavy quarks. 

The experimental data [33] at RHIC energy 
(v / ^vjv=200 GeV) shows non-zero elliptic flow in- 
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FIG. 2: Same as Fig. [T] for 0-20% centrality. Experi- 
mental data contain both the statistical and systematic 
errors, taken from [37|. 
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FIG. 3: Same as Fig. [TJfor 20-40% centrality. Experimen- 
tal data contain both the statistical as well as systematic 
errors, taken from 13711. 



dicating substantial interactions of the plasma par- 
ticles with charm and bottom quarks from which 
electrons are originated through the process: c(b) 
(hadronization) — > D(B)(dec&y) — > e+X. In Fig. Q] 
we compare the experimental data obtained by the 
PHENIX [33] collaborations for Au + Au minimum 
bias collisions at ^/sntv = 200 GeV with theoretical 
results obtain in the present work. We observe that 
the data can be reproduced by considering both ra- 
diative and collisional loss with c s = 1/ \J\. In this 
case v^ F first increases and reaches a maximum of 
about 8% then saturates for pr > 2 GeV. From the 
energy dissipation we have evaluated the shear vis- 
cosity (rj) to entropy (s) density ratio using the re- 
lation [34|: rj/s ~ 1.25T 3 /g, where q = (p 2 -) / L and 
dE/dx ~ asipr) HH, L is the length traversed by 
the heavy quark and pt is its transverse momentum. 
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FIG. 4: The variation of v" F with p T for LHC initial 
conditions for 0-10% centrality. 
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FIG. 5: The variation of v^ F with pr for low energy 
RHIC initial conditions for 0-10% centrality. 



The average value of ri/s ~ 0.1 — 0.2, slightly above 
the AdS/CFT bound [Hi]. In Fig. Q] we also display 
the results obtained when the expansion of the back- 
ground matter is switched off (dashed line). In this 
case values of the drag and diffusion co-efficients are 
evaluated at the initial temperature. The magnitude 
of t>2 is more when background is static, the reason 
for this can be understood from the following. The 
V2 controls the length of the major axis (= 1 + V2) 
and minor axis (= 1 — V2) and therefore, it also de- 
termine the value of the eccentricity as follows: 
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Eliminating e p from Eqs.[5]and [TU]we can write 
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Absence of expansion of the background matter is 
unable to smoothed out the prevailing momentum 
anisotropy ((p 2 ) ^ (p 2 ,)) in the system hence a larger 
value of vi is obtained for a static system as indi- 
cated by Eq. [TT] The difference between the dashed 
and solid lines in Fig. Q] also indicates the role of 
the expansion of the background vis-a-vis initial ge- 
ometric anisotropy. We would like to mention that 
except the results indicated by dashed line in Fig. [T] 
all other results are obtained with the background 
expansion within the ambit of (2+1) dimensional hy- 
drodynamics. 

The comparison of results obtained in the present 
work with the experimental data (37j obtained for 
two different centralities are depicted in Figs. [2] and 
[3J The agreement is reasonable both for — 20% and 
20 - 40% centralities. 

The v§ F (for 0-10% centrality collision) at LHC 
energy has been obtained by using the current for- 
malism is displayed in Fig. [4] The variation of v^ F 
with pt is similar to that of RHIC. 

Now we turn to collisions at low energy, ^/snat = 
39 GeV. For the initial conditions the pr distribu- 
tion of heavy quarks in pp collisions is required. At 
low collision energy rigorous QCD based calculations 
for heavy flavour production is not available. In 
the present work this is obtained from pQCD cal- 
culation [H, [39| for the processes: gg — > QQ and 
qq — > QQ. The composition of the matter formed 
at low energy collisions will be different from that 
of the matter formed at highest RHIC/LHC ener- 
gies. The net baryon number will be non-zero for 
low energy collisions. Therefore, in addition to tem- 
perature we need another thermodynamic variable, 
the baryonic chemical potential, fis to character- 
ize the matter. As a result in addition to energy 
momentum conservation equation which is used to 
study the evolution of baryon free matter we need 
to solve the baryon number conservation equation 
also. More details about the initial conditions, T 
and jiB dependence of the transport coefficients and 
the space time evolution have been outlined in [20J. 

At low ^/sNAr(=39 GeV) the net baryon density 
at mid-rapidity is non-zero and its value could be 
high depending on the value of v / snjv- Therefore, 
we need to solve the FP equation for both non-zero 
T and \ib which enters the equation through the 
drag and diffusion coefficients. The baryonic chem- 
ical potential for y / sn1v=39 GeV has been obtained 
from the parametrization of the experimental data 
on hadronic ratios as (see also |4l|). 



Ms(sNiv) = a(l + V^NN/bY 



(12) 



where a = 0.967 ± 0.032 GeV and b = 6.138 ± 0.399 
GeV. The parametrization in Eq. [12] gives the values 
of /is at the freeze-out. The corresponding value 
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at the initial condition is obtained from the baryon 
number conservation equation. 

In Fig. [5j v^ F for ^/snjv = 39 GeV has been dis- 
played. The lack of saturation in the trend of v^ F 
reflect insignificant B-meson contribution as well as 
the short life time of the plasma. 

So far we have discussed the azimuthal asym- 
metries of the non-photonic electron produced in 
nuclear collisions due to the propagation of the 
heavy quark in the partonic medium in the pre- 
hadronization era. However, the azimuthal asymme- 
tries of the D mesons in the post hadronization era 
(when both the temperature and density are lower 
than the partonic phase) should in principle be also 
taken into account. The suppression of the D mesons 
in the post hadronization era is found to be negligi- 
bly small [20] , indicating the fact that the hadronic 
medium (of pions and nucleons) is unable to drag 
the D mesons. In the same spirit we have neglected 
the azimuthal asymmetries due to the hadronic in- 
teractions. 

In summary, we have evaluated the azimuthal 
asymmetries (w|^ F ) using the Fokker-Planck equa- 
tion with the same set of inputs that reproduces the 
nuclear suppression, Raa measured experimentally 



at the highest RHIC energy. We found that the B 
meson suffer less flow than D mesons. The B-meson 
contribution reflects itself by the saturation of 
in both LHC and RHIC energy and lack of this trend 
in low energy RHIC run reflect the insignificant con- 
tribution of B-meson. The flow pattern at LHC is 
similar to that of RHIC. Some comments on the Raa 
vis-a-vis V2 flow are in order here. The Raa contains 
the ratio of px distribution of electron resulting from 
Au+Au to p+p collisions, where the numerator con- 
tains the interaction of the heavy quarks with the 
medium (QGP) and such interactions are absence in 
the denominator. Whereas for V2 both the numer- 
ator and the denominator contain the interactions 
with the medium, resulting in some sort of cancella- 
tion. Therefore, the Raa will be sensitive to the K 
factor but one should not expect the same kind of 
sensitivity in i>2- 
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